Trial of Galileon gravity by cosmological expansion and growth observations 
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Galileon gravity is a robust theoretical alternative to general relativity with a cosmological con- 
stant for explaining cosmic acceleration, with interesting properties such as having second order field 
equations and a shift symmetry. While either its predictions for the cosmic expansion or growth 
histories can approach standard ACDM, we demonstrate the incompatibility of both doing so simul- 
taneously. Already current observational constraints can severely disfavor an entire class of Galileon 
gravity models that do not couple directly to matter, ruling them out as an alternative to ACDM. 



I. INTRODUCTION 

General relativity is an excellent description of gravi- 
tation on all scales at which it has been tested, from the 
solar system to compact objects to cosmology. However 
within cosmology, general relativity requires a cosmologi- 
cal constant or some form of strongly negative pressure to 
explain observations of late time cosmic acceleration [1- 
5] . No compelling explanation exists for the magnitude of 
the cosmological constant (or scalar field potential), and 
in general such a contribution to the gravitational action 
should receive corrections from high energy physics. 

These fine tuning and naturalness issues motivate ex- 
ploration of further physics that can explain the observed 
acceleration and cosmic gravity, while being protected 
against high energy radiative corrections. One of the 
most successful such theories is Galileon gravity [6-8]. 
In the standard, cosmological constant free case this in- 
volves four terms in a Lagrangian that leads to well be- 
haved, second order equations of motion. The Galileon 
field arises as a geometric object in higher dimensions 
and acts in 4D like a shift symmetric scalar field (dark 
energy) , protected against radiative corrections in the ab- 
sence of matter. When introducing matter, a coupling is 
generically expected to arise which would break the shift 
symmetry and spoil the radiative stability of the model. 
However throughout this work we concentrate solely on 
a Galileon scalar field that does not couple explicitly to 
matter (we call such a model the 'uncoupled' Galileon). 

In this Letter we confront this robust theoretical al- 
ternative to general relativity with current observational 
data. This serves as a significant example of advancing 
application of cosmological data to probe gravity. One of 
the leading alternatives, Galileon gravity specifically has 
intriguing cosmological properties such as evolution from 
a high rcdshift matter dominated attractor to current ac- 
celeration and a future de Sitter attractor, and having a 
time varying effective Newton's constant G e ff • These cos- 
mological characteristics were investigated in detail in [9] 
(the existence of an asymptotic de Sitter state was first 
found in [10]); numerous other studies [11-18] have ex- 
amined various properties of Galileon gravity. Despite its 
attractive theoretical properties, though, does Galileon 



gravity remain viable obscrvationally - what is the pro- 
bative power of current data? 

We begin by an analytic discussion of the general de- 
pendencies of the effective dark energy equation of state 
and gravitational strength on the Galileon parameters. 
These will exhibit a tension between the cosmic ex- 
pansion and growth trends, so we proceed numerically 
with a full Markov Chain Monte Carlo scan through the 
Galileon parameter space, comparing the theoretical pre- 
dictions to the latest cosmological data sets. 



II. COSMOLOGICAL PROPERTIES OF 
GALILEONS 

The Galileon action is that of a scalar field ir non- 
lincarly and derivatively coupled to itself, and in curved 
spacctime to the Ricci tensor and its contractions. The 
action leads to field equations no higher than second or- 
der (and is hence a subset of Horndcski's theory [19]), and 
is invariant under the symmetry tt 
flat space limit, for constant c and b p 
results in four invariant combinations [7] 
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where R is the Ricci scalar, G vp the Einstein tensor, and 
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Mpi-Hg with M p i the Planck mass and Hq the 



Hubble constant. The full action is then 
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where C2-5 are arbitrary dimcnsionless constants, g is 
the determinant of the metric, and £ m is the matter La- 
grangian that contains no tt dependence in this work. 
Generalization of the coefficients to be functions of the 
field and its canonical kinetic term is possible [20-22], 
but we consider the standard Galileon case where the 
coefficients are constants. 
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First, let us explore the broad effects of the Galileon 
parameters. At high redshift, as discussed by [9], the ef- 
fective dark energy equation of state w(z) follows tracker 
trajectories given by the background equation of state, 
i.e. radiation or matter domination, and so is indepen- 
dent of the parameters. The strength of the gravita- 
tional coupling G e ff, however, deviates from Newton's 
constant Gn by an amount proportional to the dark en- 
ergy density at the time, fide- Thus a key early parameter 
is the initial dark energy density p Vt i{cn,Hi,Xi), where 
x = d(n I M p i) I d\n(l + z) is the field velocity and z the 
cosmic redshift. 

Analytically, increasing p n ,i increases the gravitational 
strength. A substantial increase in the gravitational 
strength would enhance the growth of structure, even at 
later times, enough to make it discrepant with observa- 
tions. So we expect that growth constraints would favor 
low values of pn,i-, keeping G c s ~ Gn for the matter 
dominated era. 

[More technically, since pv^ is a function of c„, this 
favors a certain region of the Galileon parameter space. 
At high redshift the C5 term typically dominates over the 
others, by factors of H 2 x ^> 1. For increasingly larger 
initial densities (and larger C5), it takes longer for the 
other c n terms to give comparable contributions. Since 
the moderate redshift (z ss 10) peak in the gravitational 
strength G c s noted in [9] occurs due to interplay and par- 
tial cancellation between the terms, then higher values of 
p n> i shift the peak to later times. Once the cancellation 
passes, the peak in the gravitational strength often gives 
way to a period around z ss 3 where G c ff ~ Gjv is re- 
stored. Finally, the growth of the dark energy density 
fraction Q^e moves G e s instead toward its late time de 
Sitter attractor behavior, which is independent of p n j. 
The basic point, however, is that increasing p v ^ tends to 
amplify the deviation from Einstein gravity, particularly 
at z « 3 - 10.] 

The opposite dependence is true, however, for the ex- 
pansion constraints. If we start with a low p Wl i, then due 
to the approximate tracking behavior of p n during mat- 
ter domination and the fact that we still need to arrive 
today at Qd e ,o ~ 0.7, one requires a more extreme evolu- 
tion in w(z) near the present. Analytically, to catch up 
the dark energy density to the present value one must 
have highly negative values of w(z) at < z < 2. Thus, 
low leads to strong spikes in w(z). This shifts the 
distance-redshift relation from the observed, near-ACDM 
behavior, both at low redshift and for the integrated dis- 
tance to the CMB last scattering surface. 

Thus one has a simple analytic picture: Galileon grav- 
ity is caught between the Scylla of high initial density 
(and the related region of c„ parameter space) pulling 
G e fi and growth unviably up, and the Charybdis of low 
initial density pulling w(z) and distances unviably down. 
(In Homer's Odyssey, Scylla was a cliff-dwelling monster 
pulling sailors up from ships and Charybdis a sea mon- 
ster sucking them down.) 

Figure 1 illustrates this tension between the growth 



history and expansion history behaviors in Galileon cos- 
mology. The question is whether there is a safe path 
between the monsters. This requires exact numerical 
computation, scanning over the Galileon parameter space 
with Markov Chain Monte Carlo techniques. 
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FIG. 1. Gravitational coupling deviation G b s/Gn and ef- 
fective equation of state w are shown for examples with 
high redshift initial conditions p^j = 10~ 6 p m ,i (dashed) and 
pn,i = 3 x 10~ 5 p m ,i (solid). Note that adjusting p n ,i to lessen 
deviations in gravity increases the deviations in equation of 
state, and vice versa. 



III. COSMOLOGICAL CONSTRAINTS 

We first emphasize that when using only distance data 
constraints, for example those arising from the cosmic 
microwave background (CMB), baryon acoustic oscilla- 
tion (BAO) and supernovae, we can find an acceptable 
fit, with a maximum likelihood comparable to ACDM 
(this was also found by [11]). Similarly, by using only 
growth data one can also find a viable parameter region 
(cf. the low initial density curve in Fig. 1). However, 
these two regions may be disjoint and the tension within 
the combined data constraints forces even the best fit to 
have a poor joint likelihood. 
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The second issue is that certain parts of parameter 
space are restricted theoretically due to ghosts or insta- 
bilities, as discussed in [9]. Indeed the best fit regions 
seem to tend to lie close to these because the best fits 
take advantage of the near cancellations between terms 
that can also lead to pathologies. We only apply the the- 
oretical criteria to the past behavior of the field, that is 
for z > 0, since we have no reason to rule out models 
based on their future (observationally untested) behav- 
ior. Imposing the restrictions at all times, including the 
future de Sitter state, would further constrain the allowed 
region, increasing the tension further. 

For each point in parameter space we solve for the 
effective dark energy equation of state ratio w(z) and the 
gravitational coupling G e g(z) using Eqs. (18)- (19) and 
(24) of [9]. To stay close to quantities best constrained 
by data, we use p W; j rather than Xi and fl ( ie,o rather than 
C2 as parameters, together with C3, C4, C5, and Hq. Wc 
adopt a theory or inflationary prior of a spatially flat 
universe. We carry out the Markov Chain Monte Carlo 
analysis of the full likelihood surface using CosmoMC [23] 
as a generic sampler. The likelihood is given by the sum 

£ = £cMB + ^SN + £bAO + £ growth ■ (6) 

We use the latest observational data to constrain the 
model. CMB data from WMAP7 is applied in the form 
of the covariance matrix for the shift parameter, acoustic 
peak multipole, and rcdshift of decoupling [3]. Since the 
Galileon model acts like the standard cosmology in the 
early universe these quantities basically measure the dis- 
tance to last scattering and the matter density. Distances 
from Type la super novae in the Union2.1 data compila- 
tion [24] constrain the expansion history at z « — 1.4. 
Distances from the baryon acoustic oscillation feature 
in the galaxy distribution, measured to 6 rcdshifts at 
z = 0.1 — 0.7 [25], probe a somewhat different cosmo- 
logical parameter combination. For growth constraints 
we use measurements of the growth rate from the Wig- 
gleZ survey at four redshifts z = 0.2 — 0.8 [26], and 
from the BOSS survey at z = 0.57 [27], in the form of 
their Eq. (18) 3x3 covariance matrix including expan- 
sion quantities, plus the Eq growth probe [28, 29]. Note 
that the main conditions needed to apply these growth 
data analyses to constrain a modified gravity model - 
that the standard z > 1000 matter transfer function and 
initial conditions are preserved, and that growth is scale 
independent over the relevant length scales - are satis- 
fied by the Galileon case. We only use data within this 
scale independent range, which is below the Hubble scale 
~ 4000Mpc and above the Vainshtcin scale ~ IMpc. A 
thorough analysis of non-linear effects, specifically the 
Vainshtein screening mechanism, is beyond the scope of 
this paper, although we expect them to be unimportant 
on scales relevant to linear perturbations. 

The results of the MCMC indicate the Galileon model 
is severely disfavored. The best fit yields a minimum 
Ax 2 = 31 with respect to the best fit ACDM model, de- 
spite the Galileon case having 4 extra fit parameters. We 



conclude that the entire parameter space of the standard 
Galileon theory is strongly disfavored. The CMB dis- 
tance to last scattering deviates by ~ 3a from the best 
fit ACDM case, and the individual lower redshift dis- 
tances and growth predictions are similarly in moderate 
disagreement with ACDM. The highest impact individ- 
ual constraint arises from the BOSS measurements at 
z = 0.57. This leverage bodes well for the impact of fu- 
ture redshift surveys on testing gravity on cosmic scales. 
The combination of all the data leads in aggregate to a 
poor fit. 

Figure 2 shows the \ 2 surface relative to the best fit 
ACDM result in the plane of the CMB shift parameter R 
and growth rate jo%{z = 0.57). These quantities serve as 
examples of cosmic expansion and growth, respectively. 
The A X 2 is large ( X ^ ali - 587.2 to Xacdm = 556.5 for 
the best fits) and the Galileon values are shifted con- 
siderably in attempting to fit the expansion and growth 
simultaneously (high yellow triangle for Galileon gravity 
vs low purple square for ACDM). 




FIG. 2. Ax 2 relative to the best fit ACDM model (pur- 
ple square) is shown for the Galileon model as a function of 
the CMB shift parameter R (an example of expansion) and 
growth rate fas(z = 0.57) (an example of growth, measured 
by BOSS [27]). Blue stars are points derived from the MCMC 
chains, outlining the rough paraboloid of the x 2 surface, with 
the yellow triangle the best fit. No values of the Galileon 
parameters provide a fit with A^ 2 < 30.7 (horizontal grid). 

The key tension between expansion and growth will 
be more fully realized with more accurate data. If as an 
example of future data we merely change the SN data im- 
plementation to employ the Union2.1 statistics-only error 
matrix, rather than the full systematics matrix, the im- 
proved distance measurements exhibit the tension much 
more clearly, leading to a minimum Ax 2 = 53. This 
demonstrates how upcoming supernova surveys will also 
deliver substantial cosmological leverage. 

The maximum likelihood values are highly stable with 
respect to variations in the prior ranges. The param- 
eters fl m ,o, h, an d p v ,i are well constrained (with best 
fits for the Galileon cosmology at 0.302, 0.714, and 
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ln(p W) j//9 m! i) = —11.09, respectively), but the coefficients 
C3, C4, C5 have strong covariance. The best fit to obser- 
vational data minimizes the deviations in growth and ex- 
pansion relative to ACDM, requiring a delicate balance 
among those Galileon coefficients. While the nominal 
best fits are respectively —2.10, —1.71, —1.77, there is a 
long, narrow region of degeneracy. (We actually also run 
extended ranges with logarithmic priors to ensure we are 
not missing a better fit. Likelihood indeed decreases for 
values of c„ with amplitudes much less than or greater 
than 1.) The degeneracy is moot, however, since the 
maximum likelihood is so poor. 

No point in the Galileon parameter space gives a rea- 
sonable fit to current data. Moreover, the best fit, poor 
though it is, is achieved by balancing on the edge of a 
precipice: the gravitational strength diverges in the very 
near future. To suppress deviations in growth and expan- 
sion simultaneously the Galileon terms are forced into a 
highly delicate, and temporary, cancellation. (Note the 
divergence of G c g may be ameliorated by effects beyond 
sub- horizon, linear perturbation theory.) Since we only 
applied our instability criterion to the past, where there 
is data, we do not rule out this model on theoretical 
grounds despite its Laplace instability (negative sound 
speed squared, c 2 < 0, for the dark energy perturba- 
tions) in the future. Figure 3 exhibits the gravitational 
strength and effective dark energy equation of state as a 
function of redshift for the best fit Galileon model. 



IV. CONCLUSIONS 

General relativity has passed all tests to date but lacks 
a clear explanation of the magnitude of the cosmological 
constant, or origin of dark energy, needed to account for 
cosmic acceleration. Two important, pressing questions 
are whether a sound alternative theory of gravity can 
explain this, and what leverage exists from current cos- 
mological data to test such theories. 

Galileon scalar fields, which have strong ties to higher 
dimensional gravity theories, can give rise to late time 
cosmic acceleration and possess well behaved, second or- 
der field equations with symmetries protecting against 
high energy physics renormalizations. 

We analytically identify, and numerically quantify, 
a tension, however, between Galileon predictions for 
the cosmic expansion history and growth history that 
severely disfavors Galileon cosmology. Confronting the 
entire class of standard, uncoupled Galileon theory with 
current observations demonstrates that the predictions 
are a worse fit than general relativity with a cosmologi- 
cal constant by A\ 2 > 30. If one wanted to abandon the 
theory prior of spatial flatness, [11] found that adding a 
free parameter for curvature improved the best fit x 2 by 
little more than one, and so would not have a significant 
effect on our conclusions. 

In this work we have focussed on the case in which the 
scalar field is not directly coupled to ordinary matter. In 
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FIG. 3. Gravitational strength Geff(z) and effective dark en- 
ergy equation of state w(z) are plotted for the Galileon model 
that best fits the current data. Different G ff superscripts cor- 
respond to the different modified Poisson equations in [9]. 

a previous paper [9] we studied the cosmological evolu- 
tion in the presence of an explicit coupling. This previ- 
ous work highlighted the existence of ghost and Laplace 
instabilities when a coupling is introduced, however an 
exhaustive scan of the parameter space for the coupled 
models remains to be undertaken (see [30] for work in 
this direction). 

It is striking and signficant that already with current 
data we can rule out an entire, theoretically viable class 
of extended gravity, one with several attractive features. 
We also established that forthcoming data will be able to 
strengthen these limits to A\ 2 > 50. More generally, the 
next generation of cosmological measurements will shed 
strong light on the distinction between modified gravity 
vs general relativity plus a physical dark energy, an ex- 
citing advance in understanding our universe. 
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